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摘要 : 为 了 探究 进化 模型 对 DNA 条 形 码 分 类 的 影响 , 本 研究 以 筋 灵 山 夜 蛾 科 44 个 种 的 标本 为 材料 ,获得 COI 基 
序列 。 使 用 邻接 法 (neighbor-joining) 、 最 大 简约 法 (maximum parsimony) 、 最 大 似 然 法 (maximum likelihood ) 以 及 
贝 叶 斯 法 (Bayesian inference ) 构 建 系统 发 育 树 , 并且 对 邻接 法 的 12 种 模型 、 最 大 似 然 法 的 7 种 模型 、 贝 叶 斯 法 的 
2 种 模型 进行 模型 成 功率 的 评估 。 结 果 表 明 ,邻接 法 的 12 种 模型 成 功率 相差 不 大 ， 较 稳定 ; 最 大 似 然 法 及 贝 叶 斯 
法 的 不 同 模型 成 功率 存在 明显 差异 , 不 稳定 ; 最 大 简约 法 不 基于 模型 , 成 功率 比较 稳定 。 邻 接 法 及 最 大 似 然 法 共 
A 6 种 相同 的 模型 , 这 6 种 模型 在 不 同 的 方法 中 成 功率 存在 差异 。 此 外 , 分 子 数 据 中 存在 单个 物种 仅 有 一 条 序列 
的 情况 , 显著 降低 了 模型 成 功率 , 表明 在 DNA 条 形 码 研究 中 , 每 个 物种 需要 有 多 个 样本 。 
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Abstract: The effects of evolution models on success rates of species identification via DNA barcoding 
were evaluated based on specimens of 44 Noctuidae species collected from Wuling Mountain, Hebei, 
northern China. Different phylogenetic tree reconstruction methods, including the neighbor-joining, the 
maximum parsimony, the maximum likelihood and the Bayesian inference, were used to reconstruct 
phylogenetic trees based on the COI gene sequences of these species, and success rates of different 
models (12 models for the neighbor-joining, 7 models for the maximum likelihood method, and 2 models 
for the Bayesian inference ) were evaluated systematically. The results showed that there were no 
significant differences in success rates of species identification among 12 different models for the neighbor- 
joining method, but significant differences in success rates were found among different evolution models 
for both maximum likelihood and Bayesian methods. As a comparison, the maximum parsimony method 
which is not based on evolution models obtained stable success rates of species identification. The 6 
evolutionary models shared by the neighbor-joining and maximum likelihood methods obtained different 
success rates between these two methods. In addition, we found that singletons in the database largely 
reduced success rates of species identification, suggesting that deep intraspecific sampling be required in 
the future studies of DNA barcoding. 
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species identification 


Ti SURE Noctuidae) ERAH 98 — AUR E, EHR 
已 知 3.5 J7 & ff (Poole, 1989) 。 近 年 来 ， 随 者 形态 
分 类 学 多 元 信息 的 整合 以 及 分 子 生 物 学 技术 手段 的 
高 速 发 展 , 目前 传统 夜 蛾 科 分 类 系统 发 生 了 很 大 变 
化 (Zahiri et al., 2010; Zahiri et al., 2012)。 但 是 由 
于 新 系统 还 存在 一 定 的 不 稳定 性 , 因此 本 文 以 传统 
的 分 类 系统 为 主 。 夜 蛾 科 昆 忠 绝 大 多 数 种 类 都 是 农 
林业 生产 中 的 重要 经 济 害 虫 ， 其 个 体 变异 丰富 、 近 
缘 种 间 相 似 性 较 高 ,经 稼 导致 分 类 鉴定 困难 ,特别 
是 其 包含 的 棉铃 虫 Helicoverpa armigera ( Hübner) ， 
1808 、 烟 青虫 Helicoverpa assulta ( Guenée) . JR t 2$ 
Mythimna separata-group , 甜菜 夜 蛾 Spodoptera exigua 
( Hübner) 、 地 老虎 类 Agrotis ipsilon-group, RAA E 
属 Ctenoplusia Dufay 等 分 布 范围 广 , 危害 严重 的 害 
Hy. 根据 外 部 形态 特征 更 易 导 致 鉴定 错误 和 混 消 。 
而 我 国 专家 对 夜 峨 科 的 研究 历史 较 短 ,因此 需要 投 
入 多 方面 的 调查 才能 有 效 控制 其 危害 (于 思 勤 ， 
1998; 申 效 诚 等 ,2007) 。 物 种 的 分 类 鉴定 是 对 其 
进行 科学 研究 的 基础 , 但 由 于 形态 学 鉴定 本 刁 的 局 
限 性 和 分 类 学 家 队伍 不 断 缩减 , 使 形态 分 类 学 发 展 
面临 巨大 的 挑战 (Dayrat，2005 )。 Hebert 等 
(2003a) 首次 分 机 了 200 个 杀 缘 关系 较 近 的 鲜 翅 目 
E HV I& JJ 658 bp 的 COI Fr Ez, 结果 表明 COI Fr EE 
能 够 10090 II He E x ioc 26 SEX FE EG RS [n] 4E, 
Hebert 等 (2003b ) 提出 DNA. 条 形 码 的 概念 。DNA 
REWA HEAR, 一直 致力 于 建立 标准 、 人 快捷 、 
低 成 本 的 物种 鉴定 方法 , 大 大 提高 了 物种 鉴别 的 效 
率 和 准确 性 (Hebert et al., 2003b; Smith et al., 
2005) , 并 且 在 鳞 翅 目 昆 虫 的 分 类 和 鉴定 中 得 到 了 
广泛 应 用 ( 李 青 青 等 , 2009 ) 。 

对 DNA 条 形 码 大 量 分 子 数据 包含 的 信息 进行 
有 效 的 分 析 是 研究 的 关键 ( 孝 相 林 和 张 淑 省 ,2000; 
DeSalle et al., 2005) 。 分 析 DNA 条 形 码 数据 的 方法 
大 都 集中 在 系统 发 育 分 析 的 研究 (Abdo and 
Golding, 2007) , 通过 物种 同 源 性 状 的 相似 性 比较 ， 
分 析 物 种 之 间 的 系统 发 育 关 系 , 探究 各 类 群 的 亲缘 
关系 和 进化 历史 , 最终 构建 全 部 生物 的 生命 之 树 
(Tree of Life, ToL) ( Maddison et al., 2007) 。 构 建 
系统 发 育 树 的 方法 大 体 有 以 下 几 种 , 基于 距离 的 方 
法 有 UPCMA 、 最 小 进化 法 (minimum evolution, 
ME) 和 邻接 法 (neighbor-joining, NJ) (Saitou and 
Nei，1987) 等 ; 基于 最 优 原 则 的 方法 包括 最 大 简约 


法 (maximum parsimony, MP) ( Sourdis and Nei, 
1988 ) 、 最 大 似 然 法 (maximum likelihood, ML) 
( Felsenstein, 1981; 
1997; Yang, 2002), LJ A& Ul HF Ef 3X: ( Bayesian 
inference, BI) ( Rannala and Yang, 1996; Yang and 
Rannala, 1997; Mau et al., 1999; Karol et al., 2001 ; 
Hall, 2005) 等 方法 。 

系统 发 育 树 的 构建 方法 会 预先 假定 一 个 进化 模 
型 (Penny et al., 1994) , 一般 而 言 ， 生 物化 学 性 质 
相近 的 碱 基 之 间 的 替代 频率 较 高 ; 在 DNA rp, 4 种 
碱 基 转 换 ( A 一 G, GA, COT, TC) PES EE BUT 
换 (A 一 C, AT, C>œG, GT) 以 及 它们 的 反问 取 
代 的 频率 要 高 (Fitch and Margoliash, 1967 ; Gojobori 
et al., 1982; Kocher and Wilson, 1991) , 3x JE [fg [u] 
会 影响 两 个 序列 之 间 预 期 的 分 虑 (Baxevanis and 
Ouellette, 2001) ; 由 于 核 苷 酸 序列 中 只 有 4 种 状 
态 , 在 滥 长 的 进化 过 程 中 常常 发 生 回复 和 平行 等 多 
重 置换 ,如 果 使 用 简单 的 序列 差异 估计 核 车 酸 的 蔡 
代数 , 会 使 核 音 酸 的 符 代 数 被 严重 低 佑 ,因此 必须 
使 用 模型 。 

单 参数 Jukes-Cantor (1969 ) 模 型 是 最 简单 的 核 
PRENA (Jukes and Cantor, 1969) , 该 模型 假 
设 所 有 位 点 的 核 昔 酸 蔡 换 相互 独立 ,每 一 位 点 核 音 
酸 单位 时 间 内 以 概率 入 演变 为 其 他 3 种 核 苷 酸 中 的 
一 种 。JC69 模型 假设 4 种 核 苷 酸 的 蔡 代 频率 相同 ， 
但 在 大 多 数 的 DNA. 片段 中 , 转换 的 频率 比 颠 换 高 。 
Kimura (1980) 提出 了 连续 时 间 的 二 参数 模型 
Kimura 模型 , 假设 所 有 位 点 的 核 车 酸 替换 相互 独 
XL. 且 每 一 位 点 的 核 苷 酸 转换 蔡 代 概率 为 ,斯 换 
蔡 代 率 为 B。 该 模型 的 假设 与 Jukes-Cantor 模型 假 
设 基 本 相同 , 重要 差别 是 每 个 位 点 的 核 苷 酸 蔡 代 转 
换 率 与 颠 换 率 不 相等 。 因 此 , 根据 不 同 的 核 苷 酸 蔡 
代 率 产生 了 不 同 的 碱 基 蔡 换 模型 。 随 着 研究 的 次 
A, 从 最 早 的 单 参数 模型 (JC69 ) 到 后 来 的 多 参数 模 
型 (CTR 等 ) ， 越 来 起 复杂 的 模型 被 相继 提出 ， 对 于 
其 他 的 核 车 酸 蔡 代 模 型 ， 可 参阅 Felsenstein (1981, 
1984)、Hasegawa 等 (1984)、Zharkikh(1994)、Nei 
和 Kumar(2000) 等 。 

许多 学 者 根据 影响 替代 的 各 种 因素 提出 了 不 同 
的 进化 模型 。 目 前 , 通过 不 同 模型 构建 系统 发 育 树 
对 DNA 条 形 码 分 类 影响 的 研究 还 鲜 有 报道 。 故 本 研 
完 试 图 以 筋 灵山 夜 蛾 科 44 种 DNA. 条 形 码 为 材料 , 通 
过 评估 模型 成 功率 来 对 上 述 问题 进行 初步 探讨 。 
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1 材料 与 方法 


1.1 样本 的 采集 及 鉴定 

本 研究 的 材料 鳞 翅 目 夜 蛾 类 及 毛 翅 目 样 本 采集 
于 雾 灵山 。 筋 灵山 位 于 河北 省 兴隆 县 境内 ， 东 经 
117?17' «17?35' , 北纬 40"29' ~ 40°36', 与 北京 市 
Ran HMR. FRUN RUER, IAHE, 海拔 
一 般 为 700 ~1 900 m, EMIRERRI 2 118 m, 


为 燕山 山脉 第 一 高 峰 ( 郭 文 增 等 , 2011) 。 其 恰 在 东 
北 、 华 北 和 壹 上 古 三 大 植物 区 系 的 交汇 处 , 三 大 区 系 
植物 羔 而 有 之 。 因 此 ,生态 系统 多 样 , 动 植物 种 类 
FE, 是 暖 温 市 生物 资源 宝库 ( 王 德 艺 等 , 1997) 。 

本 研究 共 采 集 夜 蛾 成 虫 样本 101 头 , 毛 翅 目 角 
H SURE ( Stenopsychidae ) 昆虫 4 头 。 取 样本 足 部 保 
存在 10090 酒精 中 , 将 样本 进行 展翅 , 制 成 标本 并 
担 照 , 取 标本 腹部 解剖 制作 外 殖 表 玻 族 以 备 物 种 鉴 
定 。 夜 蛾 分 类 鉴定 结果 详情 参照 表 1。 


R1 样本 采集 信息 及 其 COI 序列 的 GenBank 登录 序列 号 


Table 1 Specimens collecting data and GenBank accession numbers of their COI gene sequences 


标本 编号 
Specimen no. 
WLS100722. 715 
WLS100722. 413 
WLS100722. 059 


WLS100722. 121 委 夜 蛾 Athetis furvula 
WLS100722.742 委 夜 蛾 A. furvula 
WLS100722. 862 Athetis pallidipennis 


WLS100722. 082 A. pallidipennis 


WLS100722. 344 A. pallidipennis 
WLS100722. 349 


WLS100722. 081 Ed PRIR H. stygiana 
jr FAX Hypena zilla 


Moma kolthoffi 


WLS100722. 953 
WLS100722. 914 
WLS100722. 043 RR M. alpium 
WLS100722. 296 M. kolthoffi 
WLS100722. 347 
WLS100722. 534 日 雅 夜 蛾 I. japonica 

HERIR I. transversa 
贯 雅 夜 峨 transversa 


贯 雅 夜 蛾 I transversa 


WLS100722. 612 
WLS100722. 448 
WLS100722. 491 
WLS100722. 056 
WLS100722. 114 
WLS100722. 805 FÆR R. peustera 
WLS100722. 741 
WLS100722. 591 
WLS100722. 226 
WLS100722. 483 ZETA M. velutina 
ZU Ig M. velutina 
ZU Ig M. velutina 


M. monticola 


WLS100722. 353 
WLS100722. 376 
WLS100722. 229 


Species 
A ETC, Chrysorithrum flavomaculata 
暗 钝 夜 蛾 Anacronicta caliginea 


Abi HAT Plusiodonta casta 


FRIR Hypena  tristalis 


TA ET lambia transversa 


桃 剑 纹 夜 蛾 Acronicta intermedia 
EET Raphia peustera 


PEZAR Amphipyra schrenckii 
SEAHISE TE Hypocala subsatura 
2 LR Mythimna velutina 


GenBank 登录 号 
GenBank accession no. 
JX411791 
JX411792 
JX411793 
JX411794 
JX411795 
JX411796 
JX411797 
JX411798 
JX411799 
JX411800 
JX411801 
JX411802 
JX411803 
JX411804 
JX411805 
JX411806 
JX411807 
JX411808 
JX411809 
JX411810 
JX411811 
JX411812 
JX411813 
JX411814 
JX411815 
JX411816 
JX411817 
JX411818 
JX411819 
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续 表 1 Table 1 continued 


标本 编号 


Specimen no. 


WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 


860 
916 
449 
274 
838 
348 
938 
224 
570 
261 
343 
331 
453 
089 
517 
746 
248 
302 
749 
424 
291 
367 
275 
316 
371 
598 
934 
430 
341 
482 
802 
245 
373 
619 
478 
412 
505 
279 
064 
936 
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种 名 
Species 
monticola 
monticola 
monticola 
monticola 
monticola 
monticola 
monticola 
monticola 


monticola 


monticola 

钩 尾 夜 娥 Eutelia hamulatrix 
ERR E. hamulatrix 
ERR E. hamulatrix 
FERR Calyptra lata 
平 跨 索 夜 峨 C. lata 
平 跨 索 夜 峨 C. lata 
平 跨 索 夜 峨 C. lata 

C. hokkaida 

C. hokkaida 

鸟 夜 峨 Melanchra persicariae 
SE pq d Blasticorhinus unduligera 
FROR B. unduligera 
ZERKI Mocis ancilla 
XENON M. ancilla 
KARIR Cucullia kurilullia 
Narangodes argyrostrigatus 
N. argyrostrigatus 

N. argyrostrigatus 

N. argyrostrigatus 

N. argyrostrigatus 

焰 夜 蛾 Pyrrhia umbra 
KRUR P. umbra 

KRUR P. umbra 

Agrotis segetum 

H RRR Cosmia cara 
Pseudocosmia maculata 

P. maculata 

小 冠 微 夜 蛾 Lophomilia polybapta 
IRER Bryophila orthogramma 
AZUBETC X B. orthogramma 


GenBank 登录 号 


GenBank accession no. 


JX411820 
JX411821 
JX411822 
JX411823 
JX411824 
JX411825 
JX411826 
JX411827 
JX411828 
JX411829 
JX411830 
JX411831 
JX411832 
JX411833 
JX411834 
JX411835 
JX411836 
JX411837 
JX411838 
JX411839 
JX411840 
JX411841 
JX411842 
JX411843 
JX411844 
JX411845 
JX411846 
JX411847 
JX411848 
JX411849 
JX411850 
JX411851 
JX411852 
JX411853 
JX411854 
JX411855 
JX411856 
JX411857 
JX411858 
JX411859 
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续 表 1 Table 1 continued 


标本 编号 


Specimen no. 


WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 
WLS100722. 


579 
352 
588 
472 
583 
575 
015 
630 
509 
370 
462 
339 
155 
071 
504 
094 
526 
819 
124 
974 
4717 
024 
789 
975 
123 
856 
475 
943 
547 
752 
652 
469 
599 
648 
700 
768 


种 名 

Species 
兰 纹 夜 蛾 Stenoloba jankowskii 
乏 夜 峨 Niphonyx segregata 
ZARN. segregata 
ZARN. segregata 
ZARN. segregata 
FERIR Anterastria atrata 
"Hau Hadjina chinensis 
中 赫 夜 蛾 H. chinensis 
中 赫 夜 蛾 H. chinensis 
中 赫 夜 蛾 H. chinensis 


"PC H. 
"PC H. 
"PC H. 
"PC H. 
"PC H. 
"PC H. 


chinensis 
chinensis 
chinensis 
chinensis 
chinensis 


chinensis 


中 赫 夜 蛾 H. chinensis 

姬 夜 蛾 Phyllophila obliterata 
姬 夜 蛾 P. obliterata 

姬 夜 蛾 P. obliterata 

TNCS Plusilla rosalia 
ARRIR Eucarta virgo 
AER E. virgo 
RARR Polypogon gryphalis 
曲线 奴 夜 蛾 Paracolax tristalis 
曲线 奴 夜 蛾 P. tristalis 
曲线 奴 夜 蛾 P. tristalis 
BE Zanclognatha lunalis 
苹 眉 夜 峨 Pangrapta obscurata 
ZARR P. obscurata 

H TETUR I. transversa 


EROR Zanclognatha umbrosalis 


ffi Stenopsyche sp. 
ffi Stenopsyche sp. 
ffi Stenopsyche sp. 
ffi Stenopsyche sp. 


GenBank 登录 号 


GenBank accession no. 


JX411860 
JX411861 
JX411862 
JX411863 
JX411864 
JX411865 
JX411866 
JX411867 
JX411868 
JX411869 
JX411870 
JX411871 
JX411872 
JX411873 
JX411874 
JX411875 
JX411876 
JX411877 
JX411878 
JX411879 
JX411880 
JX411881 
JX411882 
JX411883 
JX411884 
JX411885 
JX411886 
JX411887 
JX411888 
JX411889 
JX411890 
JX411891 
JX411892 
JX411893 
JX411894 
JX411895 
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全 部 标本 均 于 2010 年 7 月 22 日 灯 诱 采集 自 河 北 雾 灵山 。All specimens were collected with light trap from Wuling Mountain, Hebei 
on July 22, 2010. 
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1.2 DNA 的 提取 、PCR 扩 增 及 测序 

利用 Biomed 组 织 / 细 胞 基因 组 DNA 快速 提取 
试剂 盒 提取 蛾 类 足 部 肌肉 组 织 的 基因 组 。COI 扩 增 
使 用 的 反应 体系 为 30 uL, 包括 模板 DNA 3 uL, 
ddH,O 10. 8 uL, 2 x mastermix 15 wwL， 引 物 
LCO1490 (10 umol L) 和 引物 HCO2198 (10 jumol/L) 
各 0.6 pL。 引 物 序列 为 COI 通用 引物 : 上 游 引 物 
LCO1490. 5'-GGTCAACAAATCATAAAGATATTGG- 
3'. 下 游 引 物 HCO2198. 5'-TAAACTTCAGGGTGAC 
CAAAAAATCA-3', 反应 条 件 为 94% 预 变性 5 
min; 94% 变性 30 s, 50C 退 火 30 s, 72% 延伸 30 s, 
30 个 循环 ; 最 后 7270 延伸 5 min。 获 得 的 产物 保存 
TE -20C YK o BEH 196 XU eH SE Fa DK DUI, BH 
性 结果 为 700 bp 左右 的 亮 带 。 将 阳性 PCR 产物 送 
与 北京 中 科 西 林 公 司 测序 , 得 到 101 条 夜 蛾 序列 及 
4 条 作为 外 和 群 的 毛 翅 目 序列 。 
1.3 数据 分 析 
1.3.1 构建 系统 发 育 树 : 测序 结果 为 * .abl 格式 的 
文件 , 利用 Chromas. exe 软件 查看 并 校对 序列 , RE 
套 峰 并 校正 不 可 信 的 位 点 , 使 用 BioEdit v7.0.5 软件 
中 内 置 的 CLUSTALW 多 序列 比 对 程序 进行 多 序列 
比 对 分 析 , 参数 为 默认 值 。 利 用 MEGAS. 0 软件 分 
别 采 用 MP 法 以 及 基于 Kimura 双 参 数 模型 (Kimura- 
2-parameter, K2P) ( Kimura, 1980) 采用 NJ 法 构建 
系统 发 育 树 ,系统 发 育 树 的 各 节点 文 持 率 以 序列 数 
据 集 1 000 重复 抽样 检验 的 bootstrap 值 表示 , 并 将 
大 于 50 的 bootstrap 值 标 注 于 树 图 中 。 
1.3.2 评估 模型 成 功率 : 本 文 对 雾 灵 山 夜 蛾 进行 
系统 发 育 分 析 采 用 了 4 种 不 同方 法 : NJ 法 、MP 法 、 
ML 法 和 贝 叶 斯 法 。 上 有 具体 参数 及 模型 设置 如 下 ，NJ 
法 与 MP 法 使 用 软件 PAUP * 4b10 ( Swofford, 
2003) 构建 系统 发 育 树 ; NJ 法 运行 12 种 模型 为 : 
JC69 (Jukes and Cantor, 1969; Durbin et al., 2002 ) , 
K2P( K80) (Kimura, 1980), F84 ( Felsenstein, 
1984), P (Nei and Kumar, 2000), TajNei ( Tajima 
and Nei, 1984) , TamNei ( TN93) (Tamura and Nei, 
1993; Nei and Kumar, 2000) , HKY85 ( Hasegawa et 
al., 1984; Li et al., 2000), K3P (Kimura, 1981), 
GTR ( Lanave et al., 1984; Rodriguez et al., 1990; 
Zharkikh, 1994), Total ( Swofford, 2003), Mean 
( Swofford, 2003) , LogDet ( Lake, 1994; Lockhart et 
al., 1994; Steel, 1994) ; ML 法 使 用 软件 PHYML 
v3.0 (Guindon et al., 2010 ) 构 建 系统 发 育 树 , 分 别 
运行 7 种 模型 ，Jc69，KP ( K80), F81 


( Felsenstein, 1981 ), F84, HKY85, TamNei 
(TN93 ) CTR; 贝 叶 斯 法 使 用 软件 MrBayes v3. 2 
( Ronquist and Huelsenbeck , 2003 ) 构 建 系统 发 育 树 ， 
运行 2 种 模型 : 默认 模型 ( default) 及 最 佳 模 型 
(MrModeltest) ， 默 认 模 型 即使 用 默认 模型 , 不 进行 
其 他 模型 设置 , 最 佳 模型 是 通过 Mrmodeltest 运行 ， 
基于 AlCc 原则 选择 出 来 的 最 佳 模 型 对 Mrbayes 运 
算 进 行 设置 。 

采用 严格 的 树 形 评 价 原则 计算 各 个 模型 成 功 
率 。 如 图 1(a) 所 示 , 问 询 序列 (query, Q) RAH 
种 A 序列 的 树枝 上 , 并 且 属 于 该 物种 A 的 其 他 所 
有 序列 没有 在 聚 类 到 树 图 的 其 他 树 校 上 , 那么 问 询 
成 功 ; 图 1(b) 所 示 , 问 询 序 列 Q 聚 类 到 包含 不 同 
物种 A 序列 、 物 种 B 序列 的 树枝 上 , 那么 问 询 序列 
既 不 属于 物种 A, 又 不 属于 物种 B, 问 询 不 成 功 。 


Al Al 
Fl pL 问 询 序列 Q 
A2 B 


Q assigned to A Q not assigned to any taxon 
a b 
图 1 严格 评价 模型 成 功率 的 标准 
Fig. 1 The strict evaluation standards for success rate 
A1 和 A, 分 别 代 表 物 种 A 的 不 同 序列 , B 代表 物种 B 的 序列 。Ai 
and A, stand for different sequences of species A, and B stands for the 


sequence of species B. 


对 不 同 数 据 集 采用 相应 的 抽样 次 数 ( Nue) 其 
值 为 数据 集 序 列 条 数 的 1 000 倍 , 模型 的 成 功率 
Rs 就 是 抽样 成 功 的 次 数 Ne。 除 以 样本 序列 抽 
样 总 次 数 N,,y。, 统计 成 功率 , 公式 如 下 。 

Nu =N xl1000,，R = NAN 


sample success success sample 


N 


结果 与 分 析 


本 实验 共 测 得 44 种 101 条 夜 峨 科 及 4 RER 
目的 COL 基因 序列 族 段 , 其 长 度 均 为 680 bp。 经 过 
序列 比 对 , 手工 校正 去 除 低 质 量 区 段 , EET SI 
列 长 度 为 600 bp, 所 有 的 序列 提交 至 GenBank 得 到 
相应 的 编号 ( 表 1)。 

2.1 系统 发 育 树 分 析 
本 研究 中 选择 了 4 REA H (Trichoptera) 昆虫 


10 期 


的 COI 序列 作为 外 群 ， 


采用 MEGA5.0 软件 分 别 使 


用 MP 法 以 及 NJ 法 构建 系统 发 育 树 。 构 建 的 系统 


发 育 树 见 图 2。 


夜 蛾 科 形态 分 类 学 家 根据 样本 的 外 部 形态 及 外 
生殖 融 特 征 , 将 101 头 夜 蛾 样本 分 为 37 属 44 种 。 


迟 美 妍 等 : 以 筋 灵山 夜 蛾 科 为 材料 评 佑 进化 模型 对 DNA 条 码 分 类 的 影响 
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图 2(A, B) 序 列 数据 构建 的 系统 发 育 树 可 以 看 出 ， 
NJ 树 和 MP 树 夜 蛾 的 聚 类 情况 完全 相同 , 并 且 与 从 


形态 学 


和 生理 结构 角度 进行 的 分 类 基本 一 致 ,但 也 


存在 一 些 差 异 。 分 析 结 果 显 示 每 个 属 均 聚 成 了 完全 


独立 的 单 系 , 正确 率 为 10092 ; 


; 在 种 水 平 上 , PR T H 


A WLS100722.449 Mythimna monticola B WLS100722 141 Afi segetum - 
WLS100722.261 Mythimna monticola 100 WER | 00745: 32 Athetis f fura E 
Wl OSS 3 Mytaimna monticoa WLS100722.862 Athetis pallidipennis 
WESI00/27:058 ihimna monticola 82L[ WLS100732.082 Athetis pallidipennis 
WLS100722.938 Mythimna monticola 100 WLS 100722.344 Athetis pallidipennis 
WES 100722348 Mythimna monticola 76 -WL3100722229 M fhimna monticola 
WLS 100722. 48 Mn imna monticola WLS100722.570 Mythimna monticola 
WISI, SON inna monton WES100722:449 Mythimna monticola 
WLS100722.916 Mythimna monticola 100752916 Myth imna monticol 2 
WES I00/25:338 Mythimna monticola WES 100722. 274 Wp imna monticola 

thimna velutina thimna monticol. E 
Nuper WE [0O73 e3 Narima monticola 
3 imna velutina o0 WLS100722.261 Mythimna monticola 
WES 100722. 082 Athetis pallidipennis WLS100722.860 Mythimna monticola 
0722. 344 Athetis pallidipennis WRIS iei va iis MEN 
93 00722.862 Athetis pallidipennis 84 WIS100755:3)6 je jin VENT 
S100722: 121 Athetis furvula 图 100 je imna velutina 
100722.742 Athetis furvula 61 WLS100722.226 Mythimna velutina 
X O89 Capt ra lata WLS100722.353 Mythimna velutina 
ypt 1 
WLS100722. 348 Ca alypira i lata 72—- WLS100722.802 Pyrroia umbra 
WLS100722.517 Ca Lond ie vd 100 FL WLS100722.373 ja. rhia umbra 
99 *'WLSI00722.746 Calyptra WI Id 245 Pyrr rhia umbra 
WLS100722.302 Cul tra i hokkaida WLS100722.371 Cucullia kurilullia m 
WLS100722 749 Cab boa bolta 100— WLS100722-114 Raphia peustera M 
P3100722 (94 ar ra s A = WLS100722.805 Raphia peustera 
W 04 Hadjina chinénsis | WLS100722.043 Moma alpium 
WLS100722.526 Had, ina chinensis OFL- WLS100722.296 Moma kolthoffi 
WE 的 27 29i 04 Hadjina ch chinensis AEAN 214 bo ee ee 
utelia hamulatrix 
WEISS Hedra o MUR Rem XY" 
l 1 utelia hamulatrix 
WLS100722.462 Hadiina chinensis WLS100722.579 Stenoloba jankowskii = 
WLSI 00725 400 Hadjina chinensis WLS100722.064 Bryoph ita DOANH E 
WLS100722.509 Hadiina chinensis 100—WLS100722.936 pp hita orthogramma 
Ah 108/45 5512 poat ina cainensis WS 100425- 100 Me anch "d persicartae Ea 
adj jin na c 100 alyptra hok 
Mo eo i hilg ob fiferata 图 WLS100722:749 9 Calyp pira hokkaida i 
ra lata 
ny pel Reala = SoL VISA Caf pira fata i 
alyptra lata 
WES 100723 024 Eucarta Virgo E 13 E WES100722248 Calyprra lata 
WLS 100722619 Agrotis segetum 一 S6 rL-WESI00722:305 Psn Pe od E 
WLS100722.064 Bryophita orthogramma 
WLS100722.478 Cosmia cara [d 
rb 936 Bryophita orthogramma WLS100722.583 Niphonyx segregata 
WLS100722.579 Sténoloba jankowskii 一 WLS100722-352 Niphonyx segregata 
WLS 10725 Sw 0 e 一 100 WL310072. 288 Niphoriyx segregata 
iphonyx segregata 
97 NL 100722. UM is — WLS100722.534 Jambia japonica 
mA 291 Blasticorhinus unduligera H WEE 10] IUNMO rare 
SE 0049: 2214 Meis aneia unduligera mi ods L— WLS100722:347 lambia transversa 
cis anci 
16 Mocis anci ir WLS100722.612 Jambia transversa 
WLS100722.715 Chrysorithrum flavomaculata Wi WLS100722.652 [ambia transversa 
WLS100722.591 Hypocala subsatura = WLS100722.934 Narangodes argyrostrigatus 
WLS100722.943 Zanclognatha lunalis [me] WLS100722.341 Narangodes argyrostrigatus 
wa E ep iq - o C- EST Memes imus 
WLS100722279 Lo yp ponia gophar m m WLS100722.482 Narangodes argyrostrigatus 
polybap e 
WLSI 00/44 8 lambia transversa WLS100722.056 Acronicta intermedia um 
ME IRE fm anea WU e aeli 三 
WISIN Emi inem V >47 pro — - 
Qu a mbia Pu 100 ^ WES 100722-575 Anterastría atrata m 
Dun Naran s argyrostrigatus WLS100722.943 Zanclognatha lunalis By 
Nar an ode es argyrostrigatus WLS 100722. 102 Zanclognatha umbrosalis mes 
ME 243 Narangodes argyrostrigatus WLS100722.413 Anacronicta caliginea E 
LS1007 入 204 Narangodes argyrostrigatus WLS100722.279 Lophomilia poly apta uu 
WES100/25.5 Narangodes argyrostrigatus WLS100722:975 ogon ali 
Polyp (ael. jd gryph rm 
vi 81000732388 ipto segregata VN 
WLS100722.472 Niphonyx segregata WLS100722- 123 Paracolax tristalis 
WLS100722. 354 2 Niphori x segregata WLS100722.291 Blasticorhinus unduligera E 
WLS 10072 Eutelia amu atrix WLS100722.367 Blasticorhinus unduligera 
WES 100727: 453 Eutelia hamulatrix 1007 WT $100722 275 Mocis ancilla 
WLS100722.343 Eutelia hamulatrix WLS100722.316 Mocis anci a fi 
100 
WLS100722.412 Siem ocosmia maculata B WLS100722. 391 pocala subsatura mu 
98 WLS100722.505 Pseudocosmia maculata WLS100722; ME pena tristalis m 
WLS100722.478 Cosmia cara m WLS100722. di Hypena stveiana 一 
WS100722123 Paracolax tristalis 90 WI S100722.953 DE 
WLS100722.856 Paracolax tristalis 图 95— WISI 100757 953 Hy Pl, illa rosalia — 
WLS10072 sj enh Prep : WLS100722:024 Eucarta virgo E 
ia peustera WLS100722.789 Eucart 
VISIO Raphia peustera B 100 ^ WLS100735:974 Phyllophila obliterata 
10b MU Jie E Too Le WES 100722819 Phyllophila obliterata 图 
172254 Moma oth EN I pnis oie 
S100722 802 P au d ied '» WLS100/22:155 Hadna chinensis 
SUBEST 图 100 WLS100722.462 Hadiina chinensis 
00722. 371 Cücullia kurilullia gn WLS100722.071 Ha t ie chinensis 
FER WLS100722.526 Hadjina chinensis 
WL 100722.413 Anacronicta caliginea WI S100722:094 H di hi 
VS l 745. 547 Pangrapta obscurata n . a ME chinensis 
WLS100722.752 Pangrapta obscurata WLS100722.509 Ha ane chinensis 
SIUE 741 Am, phi pyra schrenckii ex WLS100722.630 Ha dina chinensis 
5 Plustddonis ed um WLS100722.504 Ha dius chinensis 
SQUE 056 Acronicta fniermedia px WLS100722.015 Hadjina chinensis 
WLS100722.648 Steno X c 2h tenopsyche sp WLS100722.370 Hadjina chinensis 
psy psy 
WLS100722.768 Stenopsychidae Stenops che S WLS100722.715 Chrysorithrum flavomaculata W 
psy psy P. 
UEM pi quc VII gerere Meme a 
ENDO HRPHOPERCIN MDC I0 [— WES100/27 699 Stenopsychidae Stenopsyche sp. 


WLS100722.700 Stenopsychidae Stenopsyche sp. 
0.03 


图 2 基于 COI 序列 构建 的 夜 蛾 系统 发 育 树 
Fig. 2 Noctuidae phylogenetic tree based on the COI gene sequence data in this study 
A; 邻接 树 NJ tree based on COI gene; B: 最 大 简约 树 MP tree based on COI gene. 分 支 处 数值 表示 重复 1 000 次 后 的 bootstrap f ( >50% ) , 标尺 
示 遗 传 距离 。 右 边 黑色 条 带 分 别 代表 的 是 分 到 一 个 单 系 的 同 种 夜 蛾 , 红色 条 带 代 表 的 是 分 到 一 个 单 系 的 不 同 种 夜 蛾 。Numbers above branches 


are bootstrap values ( >50% ) under 1 000 replicates, and the scale bar represents the genetic distance. The species of the monophyletic clade marked by 
the black frame belong to the same species, while those marked by the red frame belong to different species. 
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JEJE ER IL japonica MÆRK M. alpium 2 种 夜 蛾 外 ， 
HR 42 个 种 聚 成 单 系 群 , 并 且 具 有 和 较 高 的 文 持 率 
(NJ 树 和 MP 树 的 支持 率 均 大 于 50)。 日 雅 夜 蛾 了 
japonica 与 贯 雅 夜 蛾 I. transversa ( Kononenko and 
Han, 2007) 为 同属 不 同 种 类 ,外 部 形态 非常 相近 ， 
仅 靠 外 部 形态 特征 区 分 困难 , 特别 是 在 外 部 形态 特 
征 受 损 时 更 是 如 此 , 可 通过 制作 外 生殖 顺 玻 族 子 以 
区 别 、 鉴 定 ; ARR M. alpium 5 M. kolthoffi 同样 
通过 分 子 系统 发 育 分 析 无 法 区 分 ,而 通过 制作 外 生 
殖 器 玻 片 得 到 鉴别 结果 (Kononenko and Han, 2007) 。 

虽然 系统 发 育 树 对 近 缘 种 鉴别 存在 一 定 的 问 
题 , 但 是 无 论 在 属 阶 元 还 是 种 阶 元 鉴别 都 有 很 高 的 
准确 率 。 
2.2 不 同 模 型 的 成 功率 

图 3(A, B) 表 示 所 有 模型 分 别 在 种 水 平 及 属 水 
平 上 的 成 功率 。 总 体 上 , 模型 的 属 水 平成 功率 要 比 
种 水 平成 功率 高 ; 在 数据 去 除 单条 序列 后 , MAR 
型 成 功率 显著 提高 。NJ 法 12 个 模型 在 不 同 的 条 件 
下 均 较 为 稳定 , 成 功率 相差 不 大 ; ML 法 7 种 模型 
不 稳定 , 模型 之 间 的 成 功率 存在 差距 ; MP 法 不 基 
于 模型 , 相对 其 他 方法 成 功率 较 高 并 且 稳定 。BI 法 
中 的 最 佳 模型 MrModeltest) ， 种 水 平 模型 成 功率 均 
相对 较 高 。 属 水 平 模型 成 功率 显示 : 在 属 阶 元 去 除 
单条 序列 的 条 件 下 ,成 功率 为 100% (95% CI; 
99.2% ~100% ) , 在 种 阶 元 去 除 单条 序列 的 条 件 
F, 成 功率 降 为 81. 78% (95% CI. 76. 28% ~ 
87.28% ) ,而 其 他 大 部 分 模型 的 成 功率 为 100% ， 
此 认为 最 佳 模型 存在 不 稳定 性 。 

图 3(C,D) 表 示 NJ 法 12 种 模型 在 种 水 平 及 属 
水 平 的 成 功率 。 图 3(C) 种 水 平 模型 成 功率 显示 12 
种 模型 之 间 差 异 不 大 。 在 全 部 夜 蛾 序列 、 去 除 属 阶 
元 单条 序列 、 去 除 种 阶 元 单条 序列 3 组 数据 条 件 
下 ,平均 差异 模型 Mean 的 成 功率 分 别 为 74. 36% 
(9596 CI. 67. 3696 ~ 81. 36% ) , 87. 68% (9596 CI; 
82. 98% ~ 92. 38% ), 10096 (95% CI; 99. 1% ~ 
100% ) ; 模型 P 的 成 功率 分 别 为 74.39% (95% CI: 
68.99% ~79.79% ) , 87. 6196 (9596 CI; 82. 9196 ~ 
92.3196) , 10096 (9596 CI; 99. 1% ~ 10096) ; 全 部 
差异 模型 Total 的 成 功率 分 别 为 74. 3696 (9596 CI: 
67.36% ~ 81. 36% ), 87. 7396 (9596 CI. 82.93% ~ 
92.53% ), 10096 (95% CI; 99. 1% ~ 10096 ); 这 3 
种 模型 的 成 功率 显示 略 高 于 其 他 9 种 模型 ,其 他 9 
种 模型 相差 不 大 。 图 3 (D) 属 水 平 的 模型 成 功率 显 
示 , 12 种 模型 成 功率 之 间 儿 乎 没有 差别 。 整 体 显 


zm NJ 法 的 属 水 平成 功率 高 于 种 水 平成 功率 。 

图 3(E, F) Jg ML 法 分 别 在 种 水 平 及 属 水 平 7 
种 模型 成 功率 。 图 中 种 水 平 以 及 属 水 平成 功率 显 
zh, ML 法 中 有 的 模型 成 功率 比较 低 或 者 存在 较 大 
的 差异 。 其 中 , 模型 CTR 的 成 功率 普遍 较 低 , 图 3 
( E) 种 水 平成 功率 显示 : 全 部 夜 峨 序列 、 去 除 属 阶 
元 单条 序列 、 去 除 种 阶 元 单条 序列 3 个 不 同 的 数据 
条 件 下 ，GTR 的 成 功率 分 别 为 72. 14% (9596 CI; 
65.3496 ~ 78. 94% ) , 73. 26% (9596 CI. 66. 96% ~ 
79.5696) , 86. 83% (9596 CI. 81. 8396 ~ 91. 8396) , 
尤其 是 后 两 者 , 明显 低 于 其 他 5 种 模型 F81, F84, 
HKY85, JC69, K80 和 TN93, 它们 的 成 功率 在 这 3 
个 数据 条 件 下 分 布 范围 分 别 为 68. 56% ~ 74. 4696 , 
85.43%~87.59% , 96.0396 ~ 100% 。 在 所 有 夜 蛾 
序列 的 数据 条 件 下 , 模型 K2P( K80 ) 的 成 功率 低 于 
其 他 6 种 模型 ,为 68.5696 (9596 CI. 68. 5096 ~ 
68. 62% ) 。 在 序列 数据 为 去 除 所 有 单条 序列 的 条 
件 下 , A F81, HKY85, J69 和 K2P(K80) 成 功率 
均 到 达 100% , 而 模型 F84 和 TN93 ( TamNei) 的 成 
功率 分 别 为 96.30% (95% CI. 93. 20% ~99. 40% )， 
97.5396 (95% CI: 94.83% ~ 10096 ) ，GTR 成 功率 最 
低 为 86.83% (9596 CI. 81. 8396 ~91.83% ) 。 

同 理 , 图 3(F) 属 水 平成 功率 显示 : 全 部 夜 蛾 序 
列 、 去 除 属 阶 元 单条 序列 、 去 除 种 阶 元 单条 序列 3 
个 不 同 的 数据 条 件 下 , 模型 GTR 的 成 功率 分 别 为 
82. 76% (95% CI. 77. 96% ~ 87. 56% ), 83. 15% 
(95% CI. 77.9596 ~ 88. 35% ), 81. 9996 (95% CI; 
76.3996 87.5996) , 明显 低 于 其 他 大 部 分 模型 ,并 
日 在 去 除 种 阶 元 单条 序列 的 条 件 下 反而 比 全 部 序 
列 、 去 除 属 阶 元 单条 序列 的 条 件 下 模型 成 功率 降 
低 ; 模型 FSI 也 存在 同样 的 情况 , 在 去 除 属 阶 元 单 
条 序列 的 分 子 数据 条 件 下 成 功率 为 97.74% (9596 
CI: 95.2496 — 10096 ) , 在 序列 数据 为 去 除 所 有 种 阶 
元 的 单条 序列 后 的 条 件 下 显示 了 与 其 他 模型 成 功率 
相反 的 趋势 : 成 功率 降低 为 81. 86% (95% CI: 
74.16% ~ 89. 56% ) ， 具 有 不 稳定 性 。 模 型 F84, 
TN93(TamNei) 也 具有 不 稳定 性 , 序列 数据 为 去 除 
所 有 单条 序列 的 条 件 下 ,其 他 模型 成 功率 达到 
100% , 而 这 两 个 模型 的 成 功率 分 别 为 96. 32% 
(95% CI. 93. 22% ~ 99. 42% ) , 97. 6296 (95% CI; 
95.0296 10096) 。 因 此 ，ML 法 模型 之 间 有 一 定 的 
差别 , 模型 GTR 成 功率 相对 很 低 , 最 不 稳定 ; 模型 
F84, TN93( TamNei), F81 成 功率 存在 不 稳定 性 ; 
模型 HKY85, JC69, K2P(K80) 成 功率 较 稳 定 。 
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图 3 不 同 进化 模型 在 不 同 条 件 下 的 成 功率 
Fig. 3 The success rates of different evolutionary models in different conditions 
A; 所 有 模型 种 水 平 All models at the species level; B: 所 有 模型 属 水 平 All models at the genus level; C; NJ 模型 种 水 平 NJ models at the species 
level; D; NJ 模型 属 水 平 NJ models at the genus level; E; ML 模型 种 水 平 ML models at the species level; F; ML 模型 属 水 平 ML models at the genus 
level. All; 全 部 夜 蛾 序列 All data; Genus: 属 阶 元 无 单条 序列 Database without genus singletons; Species: 种 阶 元 无 单条 序列 Database without 
species singletons. 通过 序列 数目 乘 以 1 000 次 的 随机 抽样 , 分 别 得 到 各 个 模型 在 不 同方 法 种 水 平和 属 水 平 的 成 功率 (95% 的 置信 区 间 ), 图 中 
单条 序列 指 单个 种 或 者 单个 属 仅 有 一 条 序列 。Under 1 000 random sampling, the success rate of different models at the species level and the genus 


level in different methods (95% CI) is obtained. In the figure, singletons mean that one species or one genus only includes one sequence. 


分 别 对 比 6 个 相同 模型 JC69，HKY85，K2P ”HKY85 较 稳 定 , 在 ML 法 、NJ 法 中 成 功率 大 致 相 
(K80), TamNei( TN93), F84 和 GTR 在 NJ 和 ML 同 ; 模型 K2P(K80), F84, TN93(TamNei) 在 ML 法 
法 两 种 方法 中 的 成 功率 。 结 果 发 现 , 模型 JC69， ”中 成 功率 低 于 NJ 法 ; 模型 GTR 在 ML 法 中 成 功率 
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远 低 于 NJ 法 。 
3 讨论 


使 用 不 同 进 化 模型 进行 系统 发 育 和 成 功率 检测 
结果 显示 ,NJ 法 中 各 个 模型 成 功率 大 致 相同 , HE 
在 不 同 的 数据 条 件 下 均 稳 定 。 因 此 NJ 法 模型 的 选 
择 对 系统 发 育 分 析 影 响 不 大 。ML 法 中 各 个 模型 的 
成 功率 之 间 存 在 差异 , 在 不 同 的 数据 条 件 下 表现 出 
不 稳定 ; ML 法 与 N 法 共有 6 个 相同 模型 ， 模 型 
JC69, HKY85 在 两 种 方法 中 成 功率 大 致 相同 (另外 
4 个 模型 在 最 大 似 然 法 中 的 成 功率 低 于 在 邻接 法 中 
的 成 功率 , 模型 GTR 尤其 明显 , 在 最 大 似 然 法 中 的 
成 功率 最 高 为 86.83% (9596 CI. 81.83%~91.83% )， 
不 仅 远 低 于 在 NJ 法 中 的 成 功率 , 与 ML 法 中 的 其 
他 模型 相 比 相差 也 很 大 ) , 并 不 是 模型 越 复 林 就 越 
准确 。 最 大 似 然 法 分 析 的 核心 在 于 替代 模型 ， 受 碱 
基 频 率 的 相等 或 不 等 、 转 换 和 其 换 速 率 的 相等 或 不 
等 、 不 变 位 点 比例 的 高 低 等 很 多 因素 的 影响 , 因此 
模型 的 正确 选择 是 ML 法 的 关键 所 在 (Steel and 
Penny, 2000; 张 原 和 陈 之 端 , 2003 ) ; MP 法 不 基于 
模型 ， 有 相对 较 高 的 成 功率 ; 贝 叶 斯 法 的 最 佳 模型 
( MrModeltest) 也 存在 不 稳定 性 。 要 根据 不 同 的 方 
法 选择 合适 的 模型 。 各 个 进化 统计 模型 构建 系统 发 
育 树 提供 了 更 加 坚实 的 统计 基础 , 但 是 不 同 的 模型 
以 及 不 同 的 方法 之 间 存 在 着 一 定 的 差异 。 因 此 一 般 
推荐 用 两 种 不 同 的 方法 构建 进化 树 ， 如 果 所 得 到 的 
进化 树 类 似 , 则 结果 较为 可 靠 。 

在 本 研究 中 , 发 现 仅 包 含 一 条 序列 的 单个 物种 
或 者 单个 属 对 进化 树 的 构建 有 很 大 的 影响 , 使 所 研 
究 模 型 的 成 功率 普遍 降低 。 成 功率 显示 , 在 所 有 没 
有 去 除 仅 有 单条 序列 物种 的 夜 蛾 分 子 数据 条 件 下 ， 
模型 在 NJ 法 、ML 法 和 贝 叶 斯 法 中 , 成 功率 普遍 较 
IR, 分 布 在 68.5696 ~ 83. 0096 E, 而 当 去 除了 所 
有 的 仅 有 一 条 序列 的 物种 的 数据 条 件 下 , 绝 大 多 数 
模型 的 成 功率 到 达 了 100% 。 因 此 , 在 系统 发 育 分 
析 时 , 需 保证 要 研究 的 各 个 物种 的 DNA 条 码 序列 
不 少 于 2 条 , 防止 数据 的 缺失 给 分 析 带 来 失误 。 数 
据 的 缺乏 和 基因 序列 的 复杂 性 是 进行 系统 发 育 分 析 
面临 的 一 大 问题 , 采样 不 全 和 选 样 不 够 准确 , 会 人 
为 增加 方法 的 假 阳 性 (Koski and Golding, 2001 ) ， 
随 着 矩阵 规模 的 增 大 , 丢失 数据 问题 变 得 也 更 加 突 
出 。 因 而 有 必要 利用 模拟 和 真实 的 数据 探索 系统 发 
育 推测 的 外 在 条 件 的 影响 , 以 及 丢失 数据 对 系统 发 


育 推 测 影 响 的 内 在 机 制 (Lemmon et al., 2009) 。 

但 是 任何 一 种 方法 都 不 可 能 完全 模拟 出 进化 的 
真正 历史 , 每 种 方法 都 有 缺点 。 因 此 应 该 将 多 种 方 
法 进行 综合 , 这 样 才 能 更 加 准确 、 快 速 地 进行 物种 
分 类 和 系统 进化 的 分 析 。 因 此 , 集 形 态 、 分 子 及 生 
态 数据 为 一 体 的 综合 学 科 分 类 法 (a multidisciplinary 
approach to taxonomy ) 将 会 使 昆虫 的 物种 鉴定 更 为 
科学 和 准确 (Wahlberg et al., 2005; Evans et al., 
2007 ; Nelson et al., 2007; Teletchea et al., 2008) , 
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